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Abstract 

The quantum dynamics of two-level systems under classical oscillator 
heat bath is mapped to the classical one of a charged particle under har- 
monic oscillator potential plus a magnetic field in a plane. The behavior of 
eigenstates and tunneling and localization are studied in detail. The broken 
symmetry condition and Langevin-like dissipative equation of motion are 
obtained. Some special dynamic features are considered. 



PACS numbers: 03.65. Ge, 03.65.Bz 



1 Introduction 



Tunneling and localization in dissipative system have attracted much attention for 
their theoretical and experimental studies have much enriched our understanding 
on the real physical world [1-12]. For macroscopic and mesoscopic quantum tunnel- 
ing phenomena two classes of processes are interested: the decay from metastable 
states and quantum coherence tunneling. The former had been extensively investi- 
gated by many authors [0], ||. The studies indicate that the tunneling probability is 
significantly suppressed at strong dissipation if we require the dissipative Langevin 
equation 

dV 

Ma + m + — = F ext (1.1) 
dq 

to be valid in semiclassical regime. The latter is a more difficult problem. A re- 
markable two- level model was given by Leggett et al. |], ^ f|, and the possibility 
of localization in the presence of dissipation was intensively studied. However, it 
is difficult to give a intuitive and uniformed description to shed light the interac- 
tion between two- level system and environment. In this paper, we shall study the 
behavior of two-level model of Leggett et al. by reducing the quantum dynam- 
ics of the two-level system to the movement of a charged particle in a harmonic 
oscillator potential under a magnetic field in a plane through a mapping. The 



method is the generalization of that discussed in ref. |T^]. In this approach the 
quasiclassical dissipative dynamics for two-level quantum evolution is mapped to 
a classical mechanical problem. Here "quasiclassical" means that we shall treat 
the environment oscillators classically and the two-level system interacted with the 
environment oscillators by quantum mechanics. 

Our main conclusion is the following: 
(1) When the acceleration of environment oscillators is small, the total energy val- 
ues(system plus environment) E c < —Hujq as increasing of the interaction between 
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the system and environment. The states are localized in the left or right of a 
symmetric double-well. 

(2) The solution is stable. 

(3) The dynamic equation for S(t) = (03) is obtained that explores the main 
physical picture of the system interacted with environment. 

2 Ground State 



The Hamiltonian for spin-boson is given by 

1 1 N v 2 a N 

H = --hAa, + - + m,4rj) + |a 3 £ c jXj , (2.1) 

j=i j j=i 

where A = 2u)q is the energy splitting between the two levels, 01, 03 are Pauli ma- 
trices and the coordinate and momentum of i-th oscillators, respectively. 
Eq (|2.1|) can be recast to the form 

1 a N f) 

H = --%Ao x + |a 3 E c ^^- )V2(aj + ^ + Hr > (2 ' 2) 

.7=1 ■? 



where H R = Y,f=i( a t a j + §Vi^ an d 



r , 1 FffUjJ 1 

[a j: aj] = <%, aj = -j={^j —Xj + % —^===pj ) , (2.3) 



The wave function obeys 



Making the transformation 



one obtains 



<9# 

ih— = H^. (2.4) 
dt v 1 



^ = e -% HRt ^, (2.5) 



ih^ = I --%Aa x + 1^ E cjie^af + e'^aj) J (2.6) 

Obviously the non-commutativity between and iff? only gives rise to the canon- 
ical transformation by = e luJjt a~j' that preserves Hr = J2j=i{bfbj + \)uJjTi- 



Keeping this canonical transformation in mind and still denoting the transformed 
coordinates by Xj we have 

ih— = | --hAa, + fas | (2-7) 

Making further transformation 

= exp{-iZ(t)°3}<f> (2-8) 



where 



we have 



or 



«W = |{E^.M* = I{^ W , (2.9) 



= (-f )(<T + e*« + ,r_e-«)0 (2.10) 



*w = (-f) e "' 1,5a - < 2 - n > 



where = ( ^ j • The normalization condition is 



a\ 2 + \P\ 2 = l. (2.12) 



Defining planar vectors 



f = Reae* x + Imae^, a = re td , 

p = Repe x + lmpe y , (3 = pe~ i<p (2.13) 



and 

1o 



fl(t) = 7(E%K = E^' ( 2 - 14 ) 
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Eq. (|2.11|) can be recast to the form: 

*f = -&-**B, (2,5) 

P = -^p+f t xB. (2.16) 

Henceforth the environment is viewed to be classical, namely any oscillator in 
the heat bath is treated as one of the driven harmonic oscillator with the driven 
force proportional to the time-dependent average value of a z through the canonical 
equation of motion: 

^ B ' + ^ B ' + ^ = ° < 2 - 17 > 

where 

S(t) = (a z ) = (*\<r 3 \V) = (0|(7 3 |0). (2.18) 
Because of eqs. ( |2.12| ), it is easy to know that 

S(t) = 2r 2 - 1 = 1 - 2p 2 , (2.19) 

(en) = r 2 6, (2.20) 

f 2 = ^ V, 'p 2 = u 2 f 2 . (2.21) 



The energy conservation reads 



r 2 +\^ 2 = r 2 > ( 2 - 22 ) 

1 2^2 1 

2' 2 0H 2 



l 'p 2 + -uolp 2 = -u 2 . (2.23) 



The total energy is 

^total = <*l#l*> 



Am,- 



= -hr 2 e + qo h(r 2 -^j:c l x l + j:^ + ^ 2 ). (2.24) 

A i=i i=i A 

Noting that the energy — Uoh(cri) is precisely the minus sign of angular momentum 
of charged particle moving in a applied magnetic field B(t) shown by eq. Q2.15Q . 



Now the two-level system obeys eq. (|2.15 ) or ( [2.16|) which looks like the dy- 
namics of charged particle experiencing harmonic force in a magnetic field per- 
pendicular to the plane in which the particle moves. This picture is similar to 



Feynman-Vernon-Hellwavth's Bloch vector stratagem JT5|. The dynamics of oscil- 
lators in heat bath is described by eq. fl2.17| ). 

Eqs. ( |2.15| ) and ( |2 . 1 7| ) allows to describe the shift of ground state of two-level 
system when the interaction between the system and environment is getting large. 

Looking at eq. (|2.15|) , we see that the circular movement orbit of charged 
particle can be verified as the ground state when Bi does not vary with time 
at r = R: 



B = Y, B 3= E 



qlc] 



2hrrinUJ 2 
3 3 J 



2R 2 ). 



Denoting R = Re r , ^V =i j = veg, we obtain from eqs. ( |2.25| ) and ( |2.15| ) 



[2R 2 - 1)(1 - 



N 



2%UJn 



o. 



i=i 



It is easy to know that 

vR = tooRVl - R? 
Hence eq. ( |2.26| ) becomes 

i? 4 



E 



rrijuj] 



0. 



The solution of eq. fl2.28| ) corresponds to the lowest-energy state of 
Denoting 

D = 1 



lQhuj 2 



E 



we find 

(1) for D < 0, no solution. 

(2) for D = 0, 



R = —,v = ±—u ,B 



0. 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



(3) for D > 0, two solutions: 

1 



Ri 



1 ± 



\ 



1/2 



(2-31) 



In the case (2) 5 = 0, \<f>) 
total energy(system plus invariant): 



_^ e ±iuj t I \ energy of system is just the 



£total = ($l^l$> = *o. 



(2.32) 



It is obvious that for the absence of interaction between the system and environ- 
ment, i.e. the magnetic field vanishes, the system is nothing but a bare two-level 
one, namely, the total energy system plus environment is equal to —huj Q , and the 
excitation state has the energy Hlvq. 
In the case (3) D > means that 

J2 N 



E 



4 — l TUjUj 



c 2 

— — > huo. 



(2.33) 



Solutions of eq. ( |2.2S| ) reads 



corresp ondingly, 



2^ 



~1F 



2, ,2 / N 



E 



i miUjf 







J 


> > 



1/2 



(2.34) 



N 



c 



1/2 



i \2 



2. .2 




(2.35) 



i=l ) 

which indicates that there exists the interaction between the system and environ- 
ment. Differing from eq. fl2.32j ) the total energy is 

JV 

"'./-/-'.; 

2" 



total 



t 1-R 2 h n 1 ^ 

-liuJa\l — V -B + ^2^"' 



R 2 



j 



a 2 N c 2 



1= 1 m jUj 



1 + 



2, ,2 IN r .2 

j 



WTui I" c 



4 



- E 



3=1 m j U j 



< -huo (2.36) 



with degeneracy. Therefore, we find the new ground state with the energy lower 
than —hujQ. When the acceleration of Xj can be neglected, the physical pic- 
ture is viewed as the following: The two-level system has the ground state with 

-^total = — ^ LUJq anc ^ ^ ne exc it & tion state with -^total = +^- a; 0; if the environment 
is completely "frozen". Whereas the ground state of the system becomes double- 
degeneracy with energy lower than —Hujq as increasing of interaction between the 
system and environment that corresponds to the non-vanishing displacements of 
the magnetic oscillators from their equilibrium points. We call the degenerated 
states localized states because in double-well system theory they represent those 
states which mainly localized in the left or right well with different parities. Later 
in this paper we shall show that the bare two-level in the system becomes unstable 
and will decay into one of the localized eigenstates with the energy lower than 
— hujo in the presence of dissipation. 

3 Stability 



Eqs.( [2.15D , ( |2.17| ), ( p.!9[ ) together with ( |2.22| ) form a set of equation determining 



both the quantum dynamics of the system and semi-classical dynamics of oscillators 
in the bath. For Bi can be neglected the eq. ( |2.25D holds. With this picture 
eq. ( 2.15| ) takes the form: 



Eq.(|J) leads to 



2 2 

¥ - r9 2 + ulr ~ ^(2r 2 - iVey = 0, (3.1) 

2h j m j u j 

2 2 

r^' + 2f^ + ||(2r 2 -l)r^^4 = 0. (3.2) 



r9(r 4 — r 2 )J = Lq = constant of motion, (3-3) 



Lq is nothing but the effective angular momentum for the charged particle moving 
in the magnetic field and 

L = -E (3.5) 

where E stands for the interacted energy of the system. With the help of eqs. (|2.30| ) 
and (|2.34j ),let us discuss the stability of solutions shown by eqs. fl2.30|) and Q2.31 
(1) For r Rj ^= we have 

9 = 



_J_ 

2r 2 - 1 



(2r 2 - l) 2 ± 



2 ' 



r 



4r 3 



u (u =f J)(2r 2 



[sr 2 - l? 



= F 1 (r). 



Since Fa(^) = 



Fi(-^= + <5r) « -4cj (cu =F J)<Jr, 



(3.6) 



(3.7) 



so that the system is stable if ujq > J, critical if ujq = J and unstable if ujq < J. 
(2) For solution given by eq. ( |2.34| ) we have r « R±, then 



^ = Z(r 2 -r 4 + ^J- 2 ) 



4J 2 (r 2 - r 4 ) 
F 2 (R± + 5r) « -8cu 2 



?o 4 



4J 2 (3r 4 -r 2 ) -cj 2 



16J 2 r 3 



Fo(r) 



(3i 



16ftV 



E 



\ j m 3^j 



5r = -2C£B ] ) 2 5r. (3.9) 



Therefore the solution R± are stable (with degeneracy). 



4 Semiclassical Equation of Motion 



To study the dynamic features we follow the standard procedure (see for example 
|l~3]1 ) by assuming N, the number of bath oscillators, is large enough so that we can 



replace the sum over j by an integration over o>j in our proceeding discussion. The 
special distribution in their path integral approach to the quantum Brownian mo- 
tion |13| will be employed. Following their discussion we shall divide the heat bath 
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effect into three parts, a renormalization part, dissipative part and random force 
one and look at how these terms affect the dynamics of tunnelling and localization. 

Let us first give a formal solution of eq. ( 2.17|) which takes the form for each i 
subscript: 

B(t) = b(t) + b (t) (4.1) 

where bo(t) = c\ cos ut + c 2 sin ojt and b{t) is a particular solution which can be 
expressed by 

b{t) = Cl {t)e lut + c 2 {t)e- iut . (4.2) 
To find a particular form of b(t) the condition 



cx(t)e wt + c 2 (t)e 



-iujt 



o 



is taken into account. It leads to the well-known solution of eq. p] 



b (t) = Bj(0)cosujt 



( smujjt\ ■ 



5,(0), 



or 



b(t) 



b(t) 



— - — sin oo j (t — t) 



dr 



rt rco 

/ drS'(r) / dujj{uj) sinu;(t — r) 
Jo Jo 



where 



Following ref. M, if the Ohmic approximation is assumed: 



then 



b(t) 



f drS(T) r due^'^ sin[(t - 
Jo Jo 



T)UJ\UJ 



uo s c +2 T{s + 2) / drS(r)(t-r) 
irh Jo 

*•(— ,— , 5 ;^2(t-r) a ) 



(4.3) 

(4.4) 
(4.5) 

(4.6) 
(4.7) 

(4.8) 



(4.9) 



where F(a, b, c; u) is hyper-geometric function. 
Let us consider a special case where s — 1, then 



(0 



h v vrA 2 + t 2 



dr 



A 



7T A 2 + (t - r) 5 



where A = c<j c 1 . When u c — > oo, A = r], eq. (|4.10| ) is reduced to 



(4.10) 



(4.11) 



This limit is nothing but the Debye distribution by making J2j ~^ Jo°° dujpiuj) []13 



[cM] s 



m w 



277CJ 2 < Vl 
> Vt 



(4.12) 



where f2 is a high frequency cut-off, r/ is the phenomenalogical friction coefficient 
in the Langevin equation eq. ( |1.1|) . 

Let Vl — > oo we have 



o(t) = — — r ■ r r"(2r — 1) 

n, ir ft 



(4.13) 



The second term is much larger than the first one for the usual environment. In 
comparison to eq. ( |2.25| ) we have 

(4-14) 

where 



C 2 



(4.15) 



With this notation the broken symmetry condition eq. (|2.33| ) is rewritten as 



2 * 

— q > huj Q 



(4.16) 



and correspondingly, 



R 



V2 



± 

B± 
^total 



2 ^ \ 



2vr 2 ^V 



1 - 



o 



=F 



4rfu; 2 „ 



2 fiVgg 
fcV 4tt 2 



2 2irhuJo 
10 



+ 



< —hujQ. 



(4.17) 

(4.18) 
(4.19) 



Obviously, eq. (|4.19| ) tells that the larger the i]Q is, the lower the energy of degen- 



erated ground state is. 

By virtue of eqs. ( |4.1| ), ( |4.4| ) and ( |4.5| ), using the continuous Debye frequency 
distribution and following we obtain 

rod 

B{t) = / p D (uj)B{uj,t)duj 
Jo 

= B s (t) + F(t) (4.20) 



where 



F(t) = B i (0) cos + sin u{t (4.21) 



is the Langevin force and carries all the characteristics of random force of classical 
Brownian motion if we assume the thermal probability distribution of heat bath 
as given by eq. ( |4.12j ). In fact, it is easy to calculate the average: 

/+oo /"+0O ui ■ 

-oo d±j Loo dx r-- ex P{-^f(i 2 3+"M)} ( 4 - 22 ) 

and find random force: 

(F(t)} = 0, 

(F(t)F(t')) = 2kT^ Sinn{t ~ t,) ^ 2kT^6(t - t') (4.23) 

where F(t) corresponds to b (t) in eq. ( [4.1| ). 

With the above knowledge we come to establish the dynamic equation. We have 
already known that the angular momentum is related to the energy of two-level 
system based on eq. ( |3.5| ). In the existence of the interaction with environment the 
energy is altered by the magnetic field, namely, the renormalized angular momen- 
tum Lr times (— %) can be viewed as the renormalized energy er of the two-level 
system. The Lr is given by 

L R (t)=r 2 6 + — ^S 2 (t) and e R = -hL R . (4.24) 
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Taking the time-derivative of Lr (or er) we get a meaningful relation similar to 
the classical Langevin dissipative dynamics: 



The first term on the right hand of eq. ( 4.25Q is the dissipative term which always 



make Lr increase, or in other words, decreased Er. The second term represents 
the work done by the Langevin force. In comparison to eq. S(t) plays the 

role similar to q(t) so we can expect a Langevin-like dynamic equation for S(t) on 
the basis of eqs. ( |2.15| ) and ( [2 . 1 7| ) , that can be written in the form: 



r -r9 2 + tolr = -r9B(t), (4.26) 
rO + 2r9 = rB{t). (4.27) 

The straight calculation gives 

^S{t) = -^ S(t) + AaL R (t)S(t) - \o?S\t) - 4 (^)L B (t)|s(t) 

where a = in Comparison with eq. (|2.1|) , the physical meaning of the RHS 

of eq. Q4.28 ) can be explained as follows: 



The first three terms are the driven forces, the fourth and fifth terms are friction 
force although they do not always make the motion retarded and the last two terms 
are random force. 

Equations (|4.25| ) , ( |4.28| ) are the central results of our quasiclassical investigation 
of dissipative dynamics of two-level system. Although the equations are still hard 
to be solved, a qualitative consideration can help us insight into some physics of this 
quantum coherent dynamics. It is not difficult to find that the renormalized angular 



momentum Lr plays important roles here. It first appears in eq. ( f4.24|) representing 



the renormalized energy of the system, then in eq. (|4 . 2 8|) acting as a time dependent 
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factor to adjust the forces acting on the charged particle. As we have studied, for 
its own properties we can find that when the broken symmetry condition is not 
satisfied, L R reaches its minimum value —uq/2 and the maximum value ui/2. On 
the opposite case in which the broken symmetry condition is satisfied L R still 
reaches its minimum value —uq/2 at the eigenstate ( |2.30| ), but now eigenstate with 

2 

L R — ujq/2 becomes a saddle point and L R reaches its maximum value ^f[^^- + 
^Y^] at the true degenerate ground states ( |2.34| ). Because the dissipation always 
makes L R increase, the dynamics with any initial state will tend to ground state 
( |2.30|) at weak coupling situation and to one of the localized ground states at 
strong coupling situation. This is consistent with our basic knowledge on the 
broken symmetry picture. 

From classical mechanics we know that any velocity independent force / can 
be written as f — —dV/dq. For the first three driven force terms in ( f4.28| ) we can 
do the same thing. The explicative potential can be written in terms of S and 
L R (t): 

V R (S, t) = \2ul - 2aL R (t)} S 2 (t) + \a 2 S A . (4.29) 

L J o 

Here V R (S,t) is time-dependent and just this time dependence exhibits the novel 
dynamic features whether or not the two-level system processes the broken symme- 
try. When flrjq 2 /2tt < hu>o, the time dependent coefficient of S 2 is never negative 
and V R (S,t) is positive everywhere except at S = 0, its only minimum point 
V R (S, t) = 0. When Qi]q 2 l /2ir > Hujq we have a critical value TiuJqK / riq^Vt for L R (t). 
Under the circumstance of L R (t) < irhujQ /flrjq^ V R (S,t) has the similar behav- 
ior. As soon as L R (t) exceed this value, the figure of V R (S,t) will be changed 
to double well form. At this situation S = becomes an unstable equilibrium 
point and the system at this point will decay due to external perturbation. At 
the same time two new degenerate minimum points appear and they will tend to 
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r± = ^yi± yl — 2n 2 h 2 LjQ/Q 2 ri 2 qQ with the increasing of L^it). These behaviors 
of Ln(t) agree with the former results of this paper. 

Interesting features also occur in the dissipative force terms in ( |4.28|) . The 
fifth term, contrary to the ordinary understanding of dissipative force , behaviors 
as an advance force and the forth term plays either a retardative or advance role 
depending on the positive or negative sign of L#(t)'s value. In other words, the 
dissipative force may not only decrease tunneling dynamics but also increase it 
depending on the initial condition of the system. 

In applying above results, we would like to discuss some special dissipative 
dynamic processes in strong coupling limit. Here we only give qualitative consid- 
eration. Detail calculations and comparison with other approaches will appear in 
elsewhere. 

First we consider the most interesting case with a fully localized initial wave 
function, saying SV(0) = 1. Then from the definition of we have L^/(0) = 
Vtriq 2 , / RitTi and the radial potential felt by the charged particle likes double well 
form with the minimum points very near ±1. So we can see that the motion can 
not leave far away from the edge of the unit circle. Since the velocity dS/dt is very 
small, the dissipative force terms in ( |4.28| ) almost do not affect the motion during 
any short time period. The particle moves like a damped oscillator and losses its 
kinetic energy with the change of L^(t). As L^j(t) reaches its maximum value, 
the particle will sit at the minimum point of the potential. 

Second let us discuss the case where the eigenstate with R = v = B = 
is taken as the initial state. In this case we can show that Su(0) = 0, and 
Lrji(0) = loq/2. Now the particle stands on the metastable point and will go out 
of it because of thermal fluctuation. 

At last suppose the eigenstate with R = v = — .8 = is taken as the 
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initial state. As discussed above we have Sm(0) = 0, and Lrju(0) = —ujq/2. 
The radial potential in this case has only one minimum and the particle is at the 
equilibrium point initially. Because both the two dissipative forces are advance 
ones now, so the range of motion will become larger and larger in the early period 
of dynamics. Therefore the particle will escape this potential well and the shape 
of potential will be changed into the double well form 

In conclusion we have studied the quasiclassical dissipative dynamics of a two- 
level system and compared it with the ordinary Langevin description. The further 
investigation along this direction and fully quantum mechanical treatment are 
deserved. 
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